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Abstract 

This paper investigates simultaneous wireless information and power transfer (SWIPT) in energy 
harvesting (EH) relay systems. Unlike existing SWIPT schemes requiring the channel state information 
(CSI) for coherent information delivery, we propose a noncoherent SWIPT framework for decode- 
and-forward (DF) relay systems bypassing the need for CSI and consequently saving energy in the 
network. The proposed SWIPT framework embraces power-splitting noncoherent DF (PS-NcDF) and 
time-switching noncoherent DF (TS-NcDF) in a unified form, and supports arbitrary M- ary nonco¬ 
herent frequency-shift keying (FSK) and differential phase-shift keying (DPSK). Exact (noncoherent) 
maximum-likelihood detectors (MFDs) for PS-NcDF and TS-NcDF are developed in a unified form, 
which involves integral evaluations yet serves as the optimum performance benchmark for noncoherent 
SWIPT. To reduce the computational cost of the exact MFDs, we also propose closed-form approxi¬ 
mate MFDs achieving near-optimum performance, thus serving as a practical solution for noncoherent 
SWIPT. Numerical results demonstrate a performance tradeoff between the first and second hops through 
the adjustment of time switching or power splitting parameters, whose optimal values minimizing the 
symbol-error rate (SER) are strictly between 0 and 1. We demonstrate that M-FSK results in a significant 
energy saving over M-DPSK for M > 8; thus M-FSK may be more suitable for EH relay systems. 
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I. Introduction 

The wireless information transfer (WIT) and wireless power transfer (WPT) technologies share 
the common origin of the electromagnetic waves, which is a direct consequence of the fact that 
radio frequency (RF) signals carry energy as well as information at the same time. Nevertheless, 
prior research activities and industrial developments on WIT have been carried out independently 
from those on WPT, due to the fact that WIT and WPT operate with rather different receiver 
power sensitivities, i.e., —60 dBm for information receivers and —10 dBm for energy receivers 
Ifn , [f2l . As a result, WIT and WPT have distinct receiver architectures and design objectives. In 
particular, WIT makes use of the active mixer for maximal information delivery rate 0 , whereas 
WPT utilizes the passive rectenna for maximal energy transmission efficiency 0. It is not until 
recently that the dual use of RF signals for simultaneous wireless information and power transfer 
(SWIPT) 0 , 0 has been studied, giving rise to an appealing solution for energy-constrained 
wireless systems such as cooperative relay networks. 

A. Related Works 

Recent advances on SWIPT have demonstrated conflicting objectives in designing WIT for 
maximal information delivery rate and WPT for maximal energy transmission efficiency, as 
witnessed by the fundamental rate-energy tradeoff 0. Varshney first introduced the notion 
of the rate-energy region for SWIPT in additive white Gaussian noise (AWGN) channels 0 , 
which was then generalized to frequency-selective fading channels 0. These pioneering studies, 
however, were based on an ideal receiver circuit which simultaneously decodes information and 
harvests energy from the same signal without any loss 0 , 0 . It has been recognized that this 
ideal receiver architecture 0, 0 may be difficult to implement given the current state-of-the-art 
of electronic circuits. A practical receiver architecture incorporating both energy harvesting (EH) 
and information decoding (ID) circuits into a single terminal was proposed in 0. Based on this 
practical architecture, the wireless terminal may operate either in the time switching mode 0, in 
which the receiver perform EH and ID in a time-division fashion, or in the power splitting mode 
ll9l . in which the receiver splits the received signal into two streams for separate EH and ID at 
the same time. The corresponding rate-energy tradeoff was studied for single-input single-output 
(SISO) 0, 0 and multiple-input multiple-output (MIMO) channels 0. The power allocation, 
beamforming design, and optimization schemes were studied to achieve various tradeoffs between 
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EH and ID. The throughput maximization for a multi-user system was addressed in [ITOl . The 
optimum power allocation maximizing the information rate while guaranteeing that the harvested 
energy is above a threshold was studied for broadband systems where perfect instantaneous 
channel state information (CSI) was available at the transmitter IflTTl . Moreover, beamforming 
scheme maximizing the harvested energy subject to information rate constraint was studied under 
imperfect CSI in ffl2ll . Similarly, beamforming design maximizing the information rate subject to 
EH constraint was addressed in lH3l . The precoder maximizing the energy efficiency in MIMO 
wiretap channel subject to secrecy rate and power constraints was studied in llT4l . These works 
|f2|, If5l-[fl4| represent typical applications of SWIPT in point-to-point (P2P) communications. 

Recently, the SWIPT has also found its application in energy-constrained wireless relay 
networks. In EH relay systems, the source and destination typically have dedicated power supplies 
from battery or power grid, while the relays are powered by the RF signals radiated by the source. 
The throughput performance of an amplify-and-forward (AF) relay network employing EH relays 
was studied for both delay-limited and delay-tolerant applications in |fT5ll . The benefits of EH 
relays and the differences between EH relay networks and conventional relay networks were 
investigated in lfl6l . In [fTTIl . a joint power allocation and relay selection scheme for AF relay 
networks with EH ability was studied. A greedy switching strategy between EH and data relaying 
was investigated for a dual-hop AF network, and the outage probability was analyzed in lfT8ll . 
Besides AF relay systems, SWIPT has also been studied for decode-and-forward (DF) relay 
systems. Specifically, a DF network consisting of multiple source-destination pairs and one EH 
relay was studied in llT9il . where the outage probability was studied taking into account the spatial 
randomness of the source-destination pairs. For the same system, the diversity performance and 
power allocation were studied in ||2Q|| . Similarly, a DF EH network consisting of one or multiple 
sources, multiple relays, and one destination was studied in [|2T|. and a more general DF system 
consisting of a multiple source-destination pairs and multiple EH relays in the presence of random 
interferences was studied in lt22l . where the small-scale fading and the large-scale path-loss were 
taken into account in the outage probability analysis. 

B. Motivation 

Due to the far-held isotropic RF radiation, the energy transfer efficiency of WPT (e.g., 
microwave power transfer) via RF signals is very low. Beamforming techniques (based on large 
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antenna arrays) can substantially boost the energy efficiency by forming focused energy beams 
using the instantaneous CSI. The existing works f5l- lf22l on SWIPT have either explicitly or 
implicitly assumed coherent WIT which requires the instantaneous CSI for information decoding. 
For EH relay systems featuring coherent WIT, the source needs to periodically send training 
symbols for channel estimation, which results in an increased network signaling overhead. 
Moreover, sophisticated channel estimation algorithm needs to be implemented in a distributed 
manner among spatially separated wireless terminals, which inevitably increases the hardware 
complexity and results in extra energy consumption. The consumed energy due to channel 
estimation and training becomes a serious issue especially for EH relay systems, as the consumed 
power reduces the net harvested energy for future data relaying. Therefore, coherent SWIPT 
requiring channel estimation might be a costly approach for energy-constrained applications. As 
a remedy, the noncoherent SWIPT which eliminates the need for training signalings and the 
channel estimation, and consequently, consumes much less energy, might be an energy-efficient 
and cost-effective solution for EH relay systems. To the best of our knowledge, however, the 
noncoherent SWIPT has not been studied in the literature, which motivated this work. 

The existing works on SWIPT have been devoted to the information-theoretic performance 
analysis in terms of the rate-energy tradeoff, ergodic capacity, outage probability, etc. However, 
insights into the SWIPT regarding the bit or symbol error rate performance associated with an 
actual detection scheme from the communication-theoretic perspective are not available. To the 
best of our knowledge, the maximum-likelihood detector (MLD) achieving noncoherent SWIPT 
has not been studied yet. In this paper, we will tackle the noncoherent MLDs which characterize 
the optimum performance benchmark for ID in EH relay systems. 

C. Contribution 

The main objective of this paper is to develop the noncoherent SWIPT framework and practical 
MLDs for EH relay networks. We consider an EH relay system consisting of a source, a desti¬ 
nation, and multiple EH relays, which can harvest energy directly from the RF signals radiated 
by the source and utilize the harvested energy to decode-and-forward the source information to 
the destination. The main contributions of the paper are summarized as follows: 

• A unified noncoherent SWIPT framework embracing both power-splitting noncoherent DF 
(PS-NcDF) and time-switching noncoherent DF (TS-NcDF) protocols is developed, which 
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supports any M-DPSK or M-FSK signalings enabling the noncoherent SWIPT. 

• Following the unified noncoherent SWIPT framework, we derive the exact MLDs for PS- 
NcDF and TS-NcDF in a unified form, for both M-FSK and M-DPSK signalings. The 
exact MLDs involve numerical integral evaluations; yet are useful for characterizing the 
optimum performance benchmark for ID in noncoherent SWIPT. 

• To reduce the computational cost of the exact MLDs, we develop closed-form approximate 
MLDs with substantially lower complexity. Numerical results demonstrate that the pro¬ 
posed approximate MLDs achieve almost identical performance to the exact MLDs, thus 
constituting a practical solution for noncoherent SWIPT. 

• Useful insights into the noncoherent SWIPT are obtained from the communication-theoretic 
perspective in terms of the symbol-error rate (SER). Specifically, we demonstrate a per¬ 
formance tradeoff between the first and second hops through the adjustment of the time 
switching coefficient or power splitting factor. Moreover, the optimal values of these system 
parameters corresponding to the minimum symbol-error rate (SER) are strictly between 0 
and 1. Finally, we demonstrate that M-FSK is an energy-efficient solution for noncoherent 
SWIPT, especially when M > 8. 

• As a byproduct of this work, we obtain generic analytical results dealing with the dis¬ 
tributions of circularly symmetric complex Gaussian (CSCG) random variables and their 
transformations. These generic results may be of general usefulness for other different 
applications. Moreover, an exact expression of the transition (error) probability is developed 
for M-DPSK. In addition, a tight closed-form approximation of the transition probability 
is developed. To the best of our knowledge, these analytical results have not been reported 
in the literature. 

D. Organization and Notation 

Section HH introduces the system model of EH relay networks. Section [TlT] develops the 
noncoherent SWIPT framework for EH relay networks. Section |TV] presents the mathematical 
fundamentals developed in this paper. Section [V] proposes the closed-form approximate as well 
as the exact MLDs. Section ED evaluates the performance of the proposed detectors and Section 
IVIII concludes the paper. 

Notation : We use A = B to denote that A is defined by B. Also, (■)*, (-) T , (-) H , || • ||, E(-), 
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and ln(-) denote the conjugate, transpose, conjugate transpose, 2-norm, expectation, and natural 
logarithm, respectively. I n denotes an nxn identity matrix, 0 denotes an all-zero column vector, 
i n denotes a column vector with 1 at its n-th entry and 0 elsewhere, and Zz denotes the phase 
of z over a 2n interval of interest. Finally, x ~ CA/"(/x, X) indicates that a? is a CSCG random 
vector with mean /v and covariance X. 

II. System Model 

Consider an EH relay system comprising a source terminal T 0 , a destination terminal T,/, 
and K relays T r , r = 1, 2, • • • , K . The source and destination are powered by dedicated energy 
sources such as battery or power grid, while the relays have no power supply and instead 
can harvest energy from the RF signals radiated by the source. Using the harvested energy, 
the relays assist the source-destination communication through the noncoherent DF relaying, 
enabling SWIPT in EH relay systems where the instantaneous CSI is not needed. 

We consider a composite fading model where the wireless channels are subject to both small- 
scale fading and large-scale path-loss. Note that the path-loss can lead to a serious deterioration 
in the power transfer efficiency and thus cannot be neglected in EH systems ll22ll . Let C tJ 
denote the path-loss associated with the channel from the transmit terminal i to the receive 
terminal j, and h %3 ~ CAf{(). 1) denotes the corresponding small-scale Rayleigh fading coefficient, 
for ij G {0d : 0r,rd}^ =1 . Note that this channel modeling implies that the actual variances of 
the fading coefficients h l3 are absorbed into the corresponding path-loss components C, v and 
thus, a general asymmetric fading model is essentially considered by choosing distinct path-loss 
components C l3 for different wireless links. We consider the noncoherent communication where 
the instantaneous CSIs, h i3 , ij G {OcZ, Or, rd}f =1 , are unknown to all terminals in the network. 
The SWIPT is accomplished through the noncoherent signalings such as M-FSK and M-DPSK. 
Specifically, the source broadcasts RF signals according to the noncoherent M-FSK or M-DPSK 
modulation. The relays utilize the harvested energy to noncoherently decode and forward the 
information to the destination, where the ultimate noncoherent detection is performed. Thus, the 
end-to-end information transfer is accomplished through pure noncoherent signalings. 

Due to the absence of the CSI, sophisticated protocols such as distributed space-time coding, 
beamforming, and relay selection which require the CSI cannot be applied. Instead, we adopt 
the time-division based protocol which coordinates multiple terminals’ transmissions in different 
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time intervals. Specifically, while the source transmits its signal, all relays stay silent and harvest 
energy from the source RF signals. After the source transmission, the relays decode the source 
information and forward it to the destination in different time intervals, allowing the destination 
to combine multiple received signals and thus enabling the diversity gain. Note that the time- 
division relaying protocol allows the developed detectors to be readily integrated into existing 
time-division multiple access (TDMA) systems, and it has been considered in numerous works 
including EH relay systems lf!31l . Ifl6ll . ifTSll — lf221 and self-powered relay networks [f23ll - [|25!l . 

For ease of the exposition, we focus on the single-antenna transmission over frequency-flat 
fading channels in this paper. However, the extension to the multi-antenna and/or frequency- 
selective fading is possible. Specifically, by applying the multi-carrier transmission such as 
orthogonal frequency division multiplexing (OFDM), the frequency-selective fading is automati¬ 
cally turned into flat fading, thus enabling the application of our scheme. Moreover, the extension 
to the multiple-input multiple-output (MIMO) configuration over each hop is also possible. The 
remaining issue is to design the noncoherent/differential space-time codes specifically suitable 
for EH relay systems, which may be considered as a future work. 

III. A Unified Noncoherent SWIPT Framework 

For EH systems, the main motivation to use noncoherent modulation schemes such as non¬ 
coherent FSK and DPSkQ coupled with appropriate noncoherent detection schemes is that 
the periodic training symbols and sophisticated channel estimation algorithms are completely 
eliminated, which results in substantial reduction of the signaling overhead and processing 
burden, and consequently reduces the energy consumption in the network. In this section, by 
modifying the power splitting and time switching receiver structures for coherent SWIPT |2l 
where the instantaneous CSI is required, we develop a unified noncoherent SWIPT framework 
embracing both PS-NcDF and TS-NcDF requiring no instantaneous CSI. 

A. PS-NcDF Relaying 

In the PS-NcDF protocol, the WPT and WIT circuits at each relay terminal function simul¬ 
taneously while the power fed to each circuit is a portion split from the front-end RF received 

'We do not consider M -QAM modulation in this paper because it is typically used for coherent communication where the 
instantaneous CSI is available. 
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Fig. 1. PS-NcDF and TS-NcDF protocols, where the solid and dashed arrow lines represent the information flow and energy 
flow, respectively. The numbers above the arrow lines indicate the portions of the power utilized for corresponding purposes. 


signal. Let T (sec) denote the total transmission time divided into K + 1 blocks, each of which 
of duration (sec) is allocated to one of the transmitter including the source and K relays, 


as illustrated in Fig. |l(a)[ Each block consists of N s symbols of duration T s (sec) each. In the 
first block, the source broadcasts its signal with power P 0 (watts). Each relay splits its received 
RF signal power into two portions: the p portion is utilized for WPT and the harvested energy 
at relay T,. is E r = ripP 0 Cor\hor\ 2 -i^p±, r = 1.• • • , K, where 0 < // < 1 is the EH efficiency 
f[2l J 2 | The remaining 1 — p portion is for WIT. Thus, the harvested power P r (watts) available 
for data relaying at T r is given by 


P r — T — 7)pPo£o r \ho r \ 2 . (1) 

K+l 

In the remaining K blocks, the relays take turns to decode and forward the source information 
to the destination using full portion of the harvested power P r in (HI). An extreme case of p = 0 
represents a pure WIT system with no EH, while p = 1 corresponds to a pure WPT system with 
no ID. In this paper, we consider a general SWIPT system featuring both WIT and WPT; thus, 
we have 0 <P< I- 

1) PS-NcDF for M-DPSK: For the M-DPSK signaling, the source message m e fA = 
{0,1, • • • , M — 1} in the l -th symbol time, l = 0,1, • • -, is conveyed via two consecutive 
information-bearing symbols s(l) and s(l — 1) according to the differential encoding s(l) = 
s(l — 1 )e- y27rr " //A/ , where s(0) = 1 is the initial reference signal. According to the power splitting, 


2 Note that the inter-relay channel (IRC) is not exploited for EH, because the energy that could be harvested at a relay from 
any other EH relays via the IRC is negligible as compared to the energy harvested directly from the source. 












































the received RF signal power at each relay is split in the proportion of p and 1 — p , where the 
fractions p and 1 — p of the power are utilized for EH and ID, respectively. Thus, the received 
signal at T r , r = 1, • • • , K, intended for ID is 

2/Or(0 = \/(l — P)PoT s £orhorS(l) + y/l — pUQ r {l) + f 0 r {l), (2) 

where T s (sec) is the symbol time, u 0r (l) ~ CM (0. : <Jq t1 ) is the AWGN introduced at the receive 
antenna at T r , and t; 0r (/) ~ CM ( 0, a^ r 2 ) is the AWGN due to the ID circuit at T r (which includes 
the RF-to-baseband conversion noise, the analog-to-digital conversion noise, circuit loss, etc) O, 
el The destination does not have a power splitter and the received signal is solely used for 
ID. Thus, the received sinal at the destination from the direct source-destination link is given by 

Vo d{l) — \J PM'sC'OdMdSi}) + Wod(0 + A) d(0> (3) 

where u 0 d(l) ~ CM(0, cr^ x ) and v 0 d(l) ~ CM(0, cr^ 2 ) are the AWGNs due to the receive 
antenna and ID circuit at T d , respectively. Each relay node assists the communication from the 
source to the destination via the DF protocol, where T r detects the source message by ML 
detection and forwards the detected message m r E Ml with full harvested power P r in (Q]) to the 
destination in a time-division manner via the M-DPSK transmission. Specifically, the transmitted 
symbols s r (l), l = 0,1, • • •, at T r are successively related according to the differential encoding 
s r {l) = Sr(l — l)e- ?27rm, "/ M , where s r (0) = 1 is the initial reference signal at T r , r = 1, • • • , K. 
The signal received at the destination from T r is given by 

Urd ("0 \/f r f s■C'rdhrd’^r'(J) T U r d[V) T (4) 

where u r d{l) ~ CM ((). a 2 rd x ) and v r( i(l) ~ CM(0, cr^dp) are AWGNs due to the receive antenna 
and ID circuit, respectively. It is assumed that all AWGNs are mutually independent. 

2) PS-NcDF for M-FSK: For noncoherent M-FSK, the source message m G Ml is transmitted 
over one of the M orthogonal carriers. The received RF signal at each receive terminal is fed 
to M orthogonal subband correlators, yielding an M x 1 baseband equivalent signal vector 
y ij = [y,j (1), • • • ,jij(M)] T , ij e {Or/, Or, rd}^L ,. The baseband equivalent signal received at 
the destination from the direct link is 


Yod VPoTsC 0 dhodlm + l Upd T Vq di 


(5) 
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where u M ~ CMii). af ul , I M ) and v od ~ CAffO. af jd 2 Im) are the AWGNs due to the receive 
antenna and the ID circuit at T d , respectively. Due to power splitting at each relay, the received 
signal intended for ID at T r . is scaled by the factor a/1 — p, and the baseband equivalent received 
signal at T r , r = 1, ■ ■ ■ , K, is given by 

yo r yj (1 //) 1 \) l' s (L{) r h() r \ n i . ] \J 1 /ni()r “h V() r , (6) 

where u 0r ~ CAf( 0, a/., i Jv/) and v 0r ~ CAf(0 , cr^, 2 / M ) are AWGNs due to the receive antenna 
and the ID circuit at T. r , respectively. Each relay decodes the source message and forwards the 
decoded message m r e M with full power P r in (O to the destination in a time-division fashion. 
The baseband equivalent signal received at the destination from T r is given by 

y rd \/P r T s Il, rd h rd lm r +l T W-rf T Vrd) (7) 

where u rd ~ CAf(0, a;. d , Im) and \ rd ~ CN(A). 2 Im) are AWGNs due to receive antenna 
and ID circuit at T,/, respectively. The AWGNs are assumed mutually independent. 


B. TS-NcDF Relaying 


In the TS-NcDF protocol, the WPT and WIT functionalities take place in a time-division 
fashion. Specifically, the total transmission time T (sec) is divided into two phases: the WPT 


phase of duration aT and the WIT phase of duration (1 — a)T, as illustrated in Fig. |l(b) 


In the WPT phase, the source radiates RF energy flow with power P 0 to the relays, where 
only the EH circuits are activated. Thus, the harvested energy at T r is E r = r]P 0 C 0r \h 0r \ 2 aT, 
r = 1 ,■■■ ,K, where 0 < r/ < 1 is the energy conversion efficiency [0. In the WIT phase, 
only the ID circuit is activated at each relay. To coordinate multiple terminals’ transmissions, 
the WIT phase is further divided into K + 1 blocks of duration each. In the first block, 

the source broadcasts information with power P 0 to all relays. The relay T r utilizes its harvested 
energy E r to decode and forward the source information to the destination in the (r + l)-th 
block, r = 1, • • • , K. Thus, the harvested power P r available for data relaying at T r is given by 

E r (A + l)7]PoCor\h 0r \ 2 a 


Pr = 


(1 ~a)T 


a 


( 8 ) 


K+l 

An extreme case of a = 0 represents a pure WIT system with no EH, and a = 1 corresponds 
to a pure WPT system with no ID. In this paper, we consider 0 < a < 1, which corresponds to 
the general SWIPT relay system featuring both WIT and WPT. 
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1) TS-NcDF for M-DPSK: In the TS-NcDF protocol, the signal received at each terminal 
during the WIT phase is solely utilized for ID, i.e., no power splitting is involved. Consequently, 
the signal model for TS-NcDF employing the M-DPSK is depicted as follows: 


Vo d(0 vPoTs^OdhodSil) + ^0 d(0 

(9a) 

Vo r(0 P b-^~s^0r^'0r < ^(0 “1“ ^0r(0 “1“ ^0r(0? 

(9b) 

2/rd(0 ^ rTsCtrdhrdSriP) ^rd(0 “1“ 

(9c) 


r = 1, • • • ,K, where ~ CJ\f(Q, cr, 2 -,) and Vij(l) ~ CJ\f{f),o^- 2 ) are AWGNs due to the 

receive antenna and ID circuit, respectively, ij G {Od, Or, rd}f =l . The source symbol s{l ) is 
related to the message m through s(l ) = s(l — l)e j27rm / M . Similarly, s r (l ) = s r (l — l)e- ?2,rmr / M . 

2) TS-NcDF for M-FSK: The signal model for TS-NcDF employing noncoherent M-FSK is 

yo d — \J PoT s £-Odhodim+l + u 0cZ + v 0 di (10a) 

YOr \/ P() P.s C Or h 0 r lm+1 + u 0r + v 0r , (10b) 

y rd 's/ Pr^s^rd^rd\m r -\-1 T Urd T V rc ;, (10c) 

for r = 1, • • • , K, where ~ a‘l j X I M ) and ~ a^ 2 I M ), ij G {0 cl, Or, rd}^ =1 . 


C. A Unified Noncoherent SWIPT Framework 


Since we consider two protocols, namely PS-NcDF and TS-NcDF, employing two different 
noncoherent modulations based on M-DPSK or M-FSK, there are in total four possible schemes 
considering different combinations, which make the derivations of the detectors rather tedious. 
For ease of development, we now unify the system models of PS-NcDF and TS-NcDF such that 
the detectors for PS-NcDF and TS-NcDF are derived in a general unified form. The key to the 
unifying process is a set of reparameterizations of the noise variances and the signal-to-noise 
ratios (SNRs) as follows: 


A 

70 d — 


PoT s Cod 

2 "7~ 2 5 

a 0d,l + a 0d,2 


2 A 2 , 2 2 A 

a kd ~ a kd, 1 + a kd,2i a 0r ~ 


(1 - p)cr^ rl + c4. 2 , PS-NcDF, 
<4,1 + <4,2, TS-NcDF, 


TO r 


A 


(l — p)P()T s Cor 

A-PU 0r,l+°0r,2 ’ 

PoTgCor 

FFTTFTZi 


PS-NcDF, 

TS-NcDF, 


A 

Trd — 


■npPoTsCorC rd 

a ld,l+ a rd,2 ’ 

a ( K +1) r/Pg T s C, o r C, rd 
A- a )Ur d}1 +^d, 2 ) ’ 


PS-NcDF, 

TS-NcDF, 


(11a) 


(lib) 


where r = 1, • • • , K and k — 0, • • • , K. 
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1) Unified PS/TS-NcDF Framework for M-DPSK: In the M-DPSK transmission, each mes¬ 
sage is conveyed via two consecutive symbols. Thus, it is convenient to represent two consecutive 
received symbols as a vector. To that end, we introduce the following vectors: 


A 

1 

T - 1 

! 

'-"O 

■«-S 

A 

s(l - 1) 

A 

S r (l ~ 1) 

Vij = 

Vij (0 

, s = 

s(l) 

1 

s r (Z) 


Using definitions in (fTIT) . we can unify PS-NcDF and TS-NcDF for M-DPSK as follows: 


VOr 0f)r (”\/ d'OrhorS ^0r) i 

(13a) 

Hod = o’od(\/dodhodS + nod), 

(13b) 

Vrd 0"rd (\/Trd| Z^Or | T n r ff , 

(13c) 


where n tJ ~ CA/"(0, J 2 ) for ij e {()d. Or, rd}^ =1 . By appropriately choosing the parameters 
according to (flTT) . the unified model (fl3l) can represent either PS-NcDF or TS-NcDF. 

2) Unified PS/TS-NcDF Framework for Noncoherent M-FSK: Based on the parameters in 
CD, the system models of PS-NcDF and TS-NcDF associated with noncoherent M-FSK can 
be expressed in a unified form as follows: 

ybr 0*0 r (\/ OOrF//+1 T H-Or) , (14a) 

yo d = %(^odim+i + n orf ), (14b) 

Yrd 0"rd ( \fdrd | ^0 r \ ^rd\m r +l T Urd) ) (14c) 

where ~ CM(f). I M ) for ij e (0d, Or, rdf and r = 1, • • • , K. For ease of readability, a list 

of symbols are summarized in Table [D 

Under the noncoherent SWIPT framework, the instantaneous CSIs hjj, ij G {0 d, Or, rd}^l 0 , are 
not available to the whole network, which makes the noncoherent detection necessary. To ensure 
the end-to-end optimality of the noncoherent communication, it is required that the destination 
as well as the relays perform the ML detection. Specifically, each relay performs its own 
(noncoherent) ML detection individually and forwards its estimated message to the destination. 
Due to potential decoding errors at the relays, the ultimate (noncoherent) ML detection at the 
destination must take into account the transition decoding error probabilities in all the relays, 
which will be treated in more detail in what follows. 
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TABLE I 

List of symbols. 


M 

Modulation alphabet size 

P r , E r 

Harvested power and energy at the r-th relay 

I\ 

Number of relays 

T s , T 

Symbol time and block time 

V 

Energy conversion efficiency 

T i 

The i-th terminal 

a 

Time switching coefficient 

hij 9 £ij 

Channel coefficient, path loss for the T,;-Tj link 

p 

Power splitting factor 

2 2 
&ij, 1» ^ij ,2 

Variances of AWGNs due to receive antenna and ID circuit 

Pa 

Source transmission power 

a iji 7 ij 

Noise variance and SNR of T,; -T-, link 


IV. Mathematical Fundamentals 

To facilitate the development of the noncoherent detection schemes, it is necessary to study 
the distributions involving combinations of CSCG and Rayleigh random variables/vectors, as 
indicated by the unified SWIPT frameworks in (fT3l) and (fT4l) . 


A. Generic Probability Analysis 


We carry out a generic probability analysis involving independent random varaibles/vectors 
as follows: 


ieo ~ CJ\f(0,0, 0 I 2 ), yo ~ CM{ 0, f2 0 / M ), 

Xi^CAf(Q,ni), i = 1,2. 


(15) 


We are interested in the following random vectors: 


Xo^X^c+Xo, 

(16) 

||> 

* 

to 

+ 

O 

(17) 


where c = [l,c] T , c is any complex constant, and 1 < p < M is any integer number. Note 
that (fT6l) and ([T71) are generic forms of y rd in (I13cl) for M-DPSK and y rd in (I14cl) for M- FSK, 
respectively. 

Lemma 1 (CSCG x Rayleigh + CSCG): Given the random variables/vectors in the PDFs 
of X 0 in (fT6l) and Yq in ( fT71 ) are, respectively, given by 


fx 0 (x) 


( \x2 — CX\ P \ 

CX P Oo(l+|cP) ) / DiD 2 n | i 2 \ \ xi + c * x 2 \ 2 \ 

71 2 fl 2 ^ 1 +|C| ^ O 0 (l + | C | 2 ) ^ ’ 

1 _imi 2 -m 2 \ Vp \*\ 

(ttDqK ° V ’ ^0 ) ’ 


(18) 


fY 0 {y) 


(19) 
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where x = [xi,x 2 ] t , y = [yi, ■ ■ ■ ,Vm] t , and If, •) is defined as 


/M) = - 




X 


-(x+A) 

I V ex/ 


dx. 


( 20 ) 


Proof: See Appendix A. ■ 

Corollary 1 (CSCG x CSCG + CSCG): Let X\ = X\X 2 c + x 0 and Y\ = X\X 2 \ p + yu. Then 
the PDF of X u f Xl (x), and the PDF of Yi, f Yl (y), are given by f Xl (x) = f Xo (x) and 

fvM = fvbiy)- 

Proof: Due to rotational invariance of the CSCG distribution lf26ll . it can be shown that 
X\X 2 follows the same distribution as Xi\X 2 \, which yields the above result. ■ 

A comprehensive study of the distribution of the product of two independent CSCG random 
variables, namely, “CSCG x CSCG”, was also presented in |[27l . However, the result in It271 was 
given in terms of modified bessel function. If one utili z es the result of [[271 . then the PDF of 
“CSCG x CSCG + CSCG” may be derived as an infinite convolution integral of the modified 
bessel function with the Gaussian PDF, which is prohibitively more complex than the result 
obtained here. Note that the finite integral of (l20l) has a well-behaved integrand involving only 
elementary functions, and it has a fast convergence speed which allows for fast evaluations at 
very high accuracy, as will be demonstrated later. 

In order to develop the MLD at the destination, the transition (error) probability of the DPSK 
detector at the relay needs to be determined. To that end, it is useful to study the distribution 
of the phase difference of two consecutively received DPSK signals. 

Lemma 2: Given the random variables/vectors in (fl5l) . we consider the 2x1 random vector 
Y 2 = [Yu Y 2 } t = X\C + where c = [1, c] T . Then the PDF of 0 = IY{Y 2 - Zc is 


fs(9) 


1 1 


2tt 1 , 7 2 |c| 2 sin 2 9 

l+7(l+|e| 2 ) 


7 c cos 6 

i a- 

arccos 

7|c| cos 6 

■ 

\ \/t 2 I c I 2 +t( 1 +I c I 2 )+ 1 / 

1 

h 

to 1 

c 

2 * 2 
z sin 

6 + 7(1 + c 

2 ) + l 


( 21 ) 


where 7 = Di/Do and 6 is defined in any 2 tt interval of interest^ 

Proof: See Appendix B. ■ 

The phase difference distribution of PSK/DPSK signals was studied for a simplistic AWGN 
channel with no fading |[28l . where the PDF was given in an integral form. In contrast, in addition 


’Mathematically, fe(0) in <211 1 is a periodic function with period 2-7 t. Thus, one may consider the domain of 6 as any 27r 
interval of interest such as [0, 27r) or [—7r,7r). 
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to the AWGN x 0 ~ CM(0, QqA), we impose Rayleigh fading where the fading coefficient is 
modeled by X 1 ~ CA/"(0, fli). For the fading scenario, the obtained PDF is given by a truly 
closed-form expression (12T1) . which is more suitable for further analysis. 

The results in Lemma [Hand Corollary Q] provide a useful characterization of the joint magni¬ 
tude and phase distribution of sophisticated Gaussian transformations, while the result in Lemma 
[2] characterizes the distribution of the input-output phase difference involving Gaussian random 
variables/vectors. These generic analytical results which have not been reported in the literature 
may be useful in a different context involving different applications. 


B. Numerical Computation of I (e, /3) in (l20l) 


To facilitate the evaluation of I(e,/3 ) in (l20l) . we now derive some alternative expressions of 
J(e, /3), which are suitable for numerical computations. 

Lemma 3 (Alternative Forms): I(e,/3 ) in (l20l) can be equivalently expressed as 


*M) 



( 22 ) 


(23) 


where K 0 (x) is the zero-order modified Bessel function of the second kind Il29l eq. (8.432.1)]. 

Proof: For / 3e < 1, (l20l) can be rewritten as J(e, (3) = y 1 [ M~ ( ' x+ ^dx—f^ M^ x+ ^dx\ . 

Similarly, I(e,/3) = ^ [ f 0 ^ \e~^ x+ ^dx-\- ^e“( x+ ^dx] for (3e > 1. For /3e = 1, we have 

I(e,/3 ) = fy M~^ x+ ^dx. These expressions have a common integral-form f“ dx 

2 2 

with a > 0. By change of variable, it can be shown that M~( x+ Mdx = J a °° 

Thus, we have f" M~^ x+ ^dx = \ / 0 °° M~( x+ Mdx = A' 0 (2a), followed by [29, eq. (3.471.9)]. 
Finally, substituting a = \Z/3/e for (3e > 1 or /3e < 1 and a = f3 for j3e = 1 yields (l23l) . ■ 

The alternative form in (l22l) indicates that the original integral of (l20l) can be seen as the 
expectation of the random transformation exp (—yyy^), where t is exponentially distributed 
with unit mean. The alternative form in (1231) has a closed-form solution only for f3e = 1. For 
f3e f 1, we are able to boil down the original integral into two terms, where the first term is 
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given in an exact closed-form while the second term is given by a sub-integral. The benefit of 
using this alternative form is that the sub-integral is taken over a smaller integral interval and 
thus can be more easily computed as compared to the original integral in (l20l) . We now exploit 
these alternative forms to derive a very useful piecewise approximation for the original integral 
with very high accuracy. 

Lemma 4: For /3e 1, J(e,/3) in ([20b can be approximated by the iV-th order expression, 

J(e,/3) = I N (e, /3) + R N , N > 1, where R N is the iV-th remainder and I N (e, /3) is given by 


/tvM) = 


re 


i 

e e 


1 

e e 


—p T 


N 


1 + £ (-1) WX £ (izM(zM 

n =2 n 

Vmi EZ~*~ J 


<?-! 


Ao2> 


Kq ( 2\/£ ] + i • 


AT 


, /3e<l,2l>l, 

/3e < 1, 4 < 1, (24) 

/3e > 1, 


where r = q = £ = £ ££'«.!, and Z + = {0,1, 2, • • • }. Also, 'k JV (-, •) is given by 

i= 1 z—0 




^Ai(a,b) = — 


N 




i— 1 


b — a a + b 
—^— z i “I-^— 


(25) 


where Wi and z t , i = 1, • • • ,N, are, respectively, the weights and zeros of the Legendre 
polynomial of order N [30[ Table 25.4], and ip(x) = 

Proof: See Appendix C. ■ 

The first piece of (f24b is a direct consequence of (l22l) followed by the Taylor series expansion 
of j-£ exp ( — y£_) around the mean E\t\ = 1. The second and third pieces of (l24l) are direct 
consequences of the Gaussian-Legendre quadrature carried out for the sub-integrals in (1231) . Note 
that the closed-form piecewise approximation of /(e,/3) by means of ijv(e,/#) in (l24l) is valid 
for any natural number N > 1, and an improved accuracy is expected by increasing the order 
N. However, from the computational cost point-of-view, it is desirable to choose a fairly small 
order N which attains the required accuracy. To that end, we set N = 2 to guarantee extremely 
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low computational complexity, and the expression of (1241) reduces to 




/3e < 1 , si. < 1, (26b) 


j3e > 1, 


(26c) 


where we have Zi = z 2 = — 3 , and w\ = W 2 = 1 for N = 2. Through extensive evaluations 


in Section [VH we will demonstrate that the piecewise approximation J 2 (e, /3) in (l26l ) is indeed 
very accurate for a wide range of e and /3 values. Considering its very low computational cost, 
we will utilize J 2 (e,/3) in the next section to facilitate the implementation of the approximate 


MLDs. 


V. MLDs for Energy Harvesting DF Networks 


Following the unified noncoherent SWIPT framework in Section |nTl we will develop unified 
noncoherent detectors in this section. Specifically, we begin by finding the exact MFDs at each 
relay, followed by the transition probability analysis. Finally, the exact MFDs as well as the 
closed-form approximate MFDs at the destination are derived for both M-DPSK and M-FSK. 

A. Unified MLDs at the Relays 

From the unified first-hop model in (U3al) and (I14al) . we see that the input-output relation is 
equivalent to the conventional P2P communication in Rayleigh fading. Thus, the MFD at the 
relay can be derived similar to that for P2P communication. Specifically, the MFDs of M-DPSK 
at T r , r = 1, • • • , K , for PS-NcDF and TS-NcDF are expressed in a unified form ETl Ch. 5.2.8] 


sf{*(0aw - ly 2 ’”"-'"}. 


(27) 


m r = arg max 

m= 0 ,--- ,M —1 


Similarly, the MFDs of M-FSK at T r for PS-NcDF and TS-NcDF are given by BTl Ch. 5.4.2] 


m r = arg 


max 

m =(),■■■ ,M —1 


y 0r (m + 1) 


(28) 
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B. Transition Probability Analysis 


Due to potential detection errors at the relay, the detected message m r at T r may be different 
from the actual source transmitted message m. In fact, a potential decision error at T r may 
transform the correct message m G M. to any incorrect message m r G A4. Thus, the transition 
probability, Pr(m r |m), namely the probability that the message m r is detected while m is actually 
transmitted serves as useful side information for the destination detection. With orthogonal 
signalings such as M-FSK, the transition from the correct message m to any one of the M — 1 
incorrect message is equiprobable OTll . Thus, the transition probability for M-FSK can be 
expressed in terms of the SER, e\, as follows: 


Pr(m r |m) = 


e r , m r = m, 




M—l ’ 


m r m, 


(29) 


where 


M—l 


£ { = 


Em) 


ra +1 


M — l 


m 


(30) 


m=1 v / 1 + m (l + 70r) ' 

For non-orthogonal signalings such as M-DPSK, however, the transition errors are not equiprob¬ 
able because each constellation point is not equidistant from all other constellation points OTIl . 
Thus, the simple relationship in (l29l) for M-FSK does not hold for M-DPSK. With non-trivial 
derivations, we derive the exact transition probability for M-DPSK as follows. 

Lemma 5 (Exact Transition Probability for M-DPSK): For the MLD of M-DPSK in (t27l) 
at relay T r , r = 1, • • • , K, the exact transition probability in Rayleigh fading is given by 
Pr(m r |m) = V\ mr - m |(7or). where 


. 2n+l . 
M y 


^n(70r) — 9 i 

111 / 2n-l . 


1 + 


7 l r sin 2 9 

1 + 2^0 r 


1 + 


70 r cos 9 


y/ 7 g r sin 2 9 + 2y 0r + 1 


arccos — 


70 r COS 9 

1 + 7dr 


d 9, (31) 


n — 0,1, ■ ■ • , M — 1. 


Proof: See Appendix D. ■ 

The exact transition probability of M-DPSK depends only on the transition phase indexed by 
\m r — m\. It is not hard to mathematically prove 'P„(7o r .) = 'Pm- nilur) for n = 0, • • • , M — 1. 
This property can be attributed to the symmetry of the PSK constellations. It suggests that 
in order to characterize M transition probabilities for M > 4, one only needs to calculate 
M/2 transition probabilities P n (7or) for n — 0, ■ ■ • , ^ — 1 via (071) . and the remaining M/2 
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transition probabilities, V n ffo r ), n = , • • • ,M, follow immediately by the symmetry P ri ( 7 o r ) = 

PM—n (7()r ) • 

Lemma 6 (Closed-Form Approx. Transition Probability): The transition probability of M-DPSK 
in Rayleigh fading can be approximated as 


Pr(m r |m) ~ < 


l-£f 

m r = m, 

£P 


2 ’ 

\m r — m\ 

o, 

otherwise, 


(32) 


where 


i 


c r 


2 ( 1 + 70 -)’ 


1.03 J l -±^L 

2cOS M 


1 - 


(^COS -§-) 


70 r 


1 + 


(i-GQS-g.) 


70 r 


M = 2, 
M > 4. 


(33) 


Proof: See Appendix E. ■ 

The approximate transition probability in (l32l) exploits the fact that the error transmissions 
from the transmitted PSK phase to its two closest neighbouring phases are most likely and 
euqi-probable. Moreover, error transitions to any farther constellations points other than the two 
closest constellation points are negligible. 


C. Unified MLDs at the Destination 

In this subsection, we develop the MLDs at the destination for both PS-NcDF and TS-NcDF 
protocols in a unified form, which exploits explicitly the transition probabilities at the relays to 
ameliorate the error propagations from the relays. 

Theorem 1 (Exact MLDs): The exact MLDs of M-DPSK at the destination for both PS-NcDF 
and TS-NcDF are given in a unified form as follows: 

27 " -Mm) + f> 

_ m r =0 


m = arg max 

m= 0 ,-" ,M— 11 l + 2y od 


r =1 


M—l 


P\m r -m\ (7or)e -/3 ' r(mr) / (2 7r J+(m r )) 


■, (34) 


where ft(m r ) ± * lyUD-VrAi-D^^ 


2 cH 


and V\ mr -m\ (I'or) is the transition probability given in OTI) . Similarly, the exact MLDs of M-FSK 
at the destination for PS-NcDF and TS-NcDF are given in a unified form as 

7o d |yod(m+l)|- 

.2 

Od 'T-— 1 


K 


m = arg max 

m=0,-,Af-l 1 + 7orf 


cr. 


]C ln 


(! -4)e 


llyrdll 2 -lyj-d( m + 1 )l 2 

a rd 
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Xl 


IrcLi 


|y rd(m + 1 )|- 


cr; 


rd 


+ 


M - 1 


M —1 

E 1 

m r =0 

rriry^m 


llyrdll 2 -|yrd( m ^ + 1 )| 2 

Hd I 


Irdi 


\jrd(m r + l)p 


a; 


rd 


•• (35) 


Proof: See Appendix F. ■ 

The exact (optimum) MLDs in (l34l) and (1351) require the evaluations of the integral I(e. if) 
in (l20l) . In addition, the transition probability involved in the exact MLD of (l34l) is given by 
a finite integral-form in (13TI) . which makes the exact MLDs complex to implement in practice. 
Nevertheless, it is very meaningful to study the exact MLDs as they serve as the optimum 
performance benchmark for SWIPT in noncoherent DF systems. With the exact MLDs, one may 
use them as a theoretical performance upper bound to evaluate the performance of suboptimum 
detectors, if developed. In the following, to make the MLDs suitable for practical systems, we 
develop closed-form approximations for the MLDs. 

Theorem 2 (Approximate MLDs): For M-DPSK, the unified-form approximate MLDs for both 
PS-NcDF and TS-NcDF are given by 


arg max < . 2 Z 0d , 

m=0,--- ,M—1 I 1+2 '1'0 d 


/do(m) + 52 ln (! - £ r) e 0r rd ,(df(m)) 

r=1 


+£ p e p T (m+1) J 2 /?+(m + 1)) 


m ~ 


arg max < . 2 Z 0d , 
m=0,---,M~l ] 1+27m 


/do(m) + 52 In (1 - ^r) e ^ {rn) h(2~i rd ,/3+(m)) 

r =1 




e -p r (m+l)j 2 ( 27rdii a+( m+1 )) +e -^ (^- 1 )/ 2 (2 7rd ,/3+(m-l)) 


M = 2, 


M > 4. 


(36) 


For M-FSK, the approximate MLDs for both PS-NcDF and TS-NcDF are unified as 


m m arg max 

m=0,— ,M —1 1 + 7od 


7o d |yod(m +1)| : 


K 


cr. 


0 d 


S hl 

r= 1 


(! -4)e 


llyrdll 2 -|yr<i( m + 1 )l 2 


X /2 ( ^rdi 


|y rd (m + 1)|- 


cr. 


rd 


J M_1 llyrdU 2 -|yrd(™r + l )| 2 

l —■ t' X ^ ‘2 


H—-— 

M- 1 ^ 


rd h Jrd, 


m r =0 

rrir^m 


|y rd{m r + 1)| 2 


cr. 


rd 


(37) 


Proof: Following the exact MLDs in Theorem Q] applying the approximation of I(e,(d ) ~ 
J 2 (e,/3) for M-FSK and M-DPSK and the transition probability approximation in Lemma |6] for 
M-DPSK immediately gives the approximate MLDs. ■ 

The ML decision metrics for the exact MLDs in Theorem Q] are given by the weighted sum 
of the decision statistics associated with the direct link and the relay branches. For each relay 
branch, the weight function is proportional to the transition probability and the decision statistic 
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is involved in the integral function J(e, /3). For the approximate MLDs in Theorem [2] the exact 
integral function I(e,/3) is replaced by its tight closed-form approximation I 2 (e,/3). In addition, 
for M-DPSK, the approximate transition probability is used in place of the exact transition 
probability, which reduces the sum of M — 1 transition terms to 2 terms, irrespective of M. 


VI. Numerical Results 


In this section, we carry out Monte Carlo simulations to gain useful design insights into the 
proposed PS-NcDF and TS-NcDF protocols. The EH efficiency at the relays is set to r] = 0.6 as 
in CD. We adopt the bounded path-loss model as in OTil in which the path-loss is strictly larger 
than one. Specifically, the path-loss from transmit terminal i to receive terminal j is modeled by 
Cij = . , where o is the path-loss exponent and D tJ is the transmitter-receiver distance for 

ij G {0d, Or, rd}^ =1 . The path-loss exponent is set to o = 2.7 which corresponds to the urban 
cellular communication environment lf32l . The source-destination distance is fixed to D od = 3 
(m), while the source-relay distance D 0r and relay-destination distance D rd can vary arbitrarily 
between 0 and 3 subject to D 0r + D rd = 3, for r = 1, • • ■ , K. For ease of simulations, the noise 
variances are set to be a‘f ] X = af ]2 = (?l/ 2 , such that each receive terminal has a total noise 
variance Oq, i.e., af j l + af j 2 = o q, where afj A and cr^ 2 , ij G {0d, Or, rd}^ =1 , are the variances 
of the receive antenna noise and the ID circuit noise, respectively. The error performance is 
parameterized by the transmitter SNR, defined as SNR = / o/cr ( j. 

We will evaluate the error performance of PS-NcDF and TS-NcDF at the same information rate 
R = — ^ — (bps), where N s is the total number of information-bearing symbols transmitted 
over T (sec). Note that N s new symbols are transmitted in the first block, which are relayed in 
the following K blocks by the relays. The symbol time T s is related to R as 



(38) 


For simplicity, we will fix N s /T = 1 such that R = log 2 M, which depends only on M. 

A. Validation of the Piecewise Approximation I 2 (e, j) in (l26l) 

We first check the accuracy of the proposed piecewise approximation J 2 (e,/3) in (l26l) . As a 
benchmark, we use the MATLAB numerical integration function “integral” to directly compute 
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Fig. 2. Piecewise closed-form approximation I 2 V, /?) vs the exact numerical integration I(e, /3) for various e and /? values. 


I(e,(3) in (l20l) at a relative error tolerance of 10 -12 , which is referred to as the “exact” result. 
We evaluate I 2 (e,/3) for ;3 e [—50,50] dB under various e values of interest ranging from 
—10 to 20 dB in Fig. [2} Since /3 spans over a very wide range of values for each e, the three 
subcases involved in (l26l) are all considered in Fig. [2] for each e value. In this case, the piecewise 
approximation of (l26l) is applied according to the conditions for the subcases. It can be observed 
that I 2 (e,/3) and J(e,/3) overlap for all possible values of e and /3, which confirms that J 2 (e, /3) 
in (l26l ) is indeed a very good approximation. 


B. Validation of the Approximate MLDs in (l36l) and (l37l) 

We validate the proposed closed-form approximate MLDs in (l36l) and (l37l) through compar¬ 
isons to the exact (integral-form) MLDs in (l34l) and (l35l) . for binary DPSK (BDPSK) and binary 
FSK (BFSK) signalings with fixed symbol rate R = 1 (bps). Specifically, the TS-NcDF protocol 
with time switching coefficient a = 0.5 in a single-relay (K = 1) EH system with D 0r = 1.5 


(m) is considered in Fig. |3(a)[ and the PS-NcDF protocol with power splitting factor p = 0.5 
for a two-relay (K = 2) EH system with D m = 1 (m) and D 02 = 2 (m) is considered in 


Fig. |3(b)[ As demonstrated in these figures, the SERs of the approximate MLDs are in excellent 
agreement with those of the exact MLDs over the whole SNR range, indicating that the proposed 
closed-form approximate MLDs are useful practical solutions for noncoherent SWIPT. 


C. Comparison of PS-NcDF and TS-NcDF 

We now compare PS-NcDF and TS-NcDF in terms of the SER at the same information rate R. 
Fig. |4] shows the SER versus the time switching coefficient a for TS-NcDF or power splitting 
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(a) (b) 

Fig. 3. SERs of MLDs for binary signalings (M = 2) with R = 1 (bps). |(a)1 TS-NcDF with a = 0.5 in a single-relay network 
( K = 1) with Dor = 1.5 (m). |(b)| PS-NcDF with p = 0.5 in a two-relay network (K = 2) with Dqi = 1 (m) and D 02 = 2 (m) 




(a) (b) 

Fig. 4. Comparison of PS-NcDF and TS-NcDF for fixed information rate in a single-relay network (K = 1) with Dor = 1.5 
(m): |(a)| M = 2, SNR = 30 dB, and R = 1 (bps); and |(b)| M = 8, SNR = 40 dB, and R = 3 (bps). 


factor p for PS-NcDF in a single-relay (K = 1) network with D 0r = 1.5 (m), where binary 


signalings (M = 2) with R — 1 (bps) and SNR = 30 dB (see Fig. 4(a)) and octonary signalings 


(M = 8) with R — 3 (bps) and SNR = 40 dB (see Fig. |4(b)| ) are considered. As shown 
in Fig. [U the SERs exhibit some global minima around a = 0.4 for TS-NcDF and around 
p = 0.8 for PS-NcDF. This is because there exists a performance tradeoff between the first 
and second hops through the adjustment of a (for TS-NcDF) or p (for PS-NcDF). Specifically, 
when a (or p) increases, more power (proportional to a or p) can be harvested and utilized 
by the second-hop data relaying, resulting in an increased SNR over the second-hop. On the 
other hand, for TS-NcDF, as a increases, less time (proportional to 1 — a) is available for the 
first-hop information transmission. Thus, the symbol time T s for TS-NcDF must decrease so as 
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(a) M = 2 


(b) M = 4 




(c) M == 8 (d) M =16 

Fig. 5. SERs of PS-NcDF and TS-NcDF employing near-optimal EH parameters with p = 0.8 and a = 0.4 for a single relay 
network (K = 1) with Do r = 1.5 (m): |(a)| M = 2; |(b)| M = 4; |(c)| M = 8; and |(d)| M = 16. 


to maintain the same target information rate R, which results in smaller received symbol energy 
and consequently decreased SNR over the first-hop. Similarly, for PS-NcDF, as p increases, less 
power (proportional to 1 — p) is split into the WIT circuit for information decoding; thus, the 
first-hop SNR decreases. Therefore, adjusting the values of a (or p) results in different impacts 
on the SNRs over the two hops, and consequently, a tradeoff exists between the first- and second- 
hop SNRs through the adjustment of a (or p). This can be also explained mathematically. By 
substituting (l38l) into (II lbl) . we see that 7 0r is monotonically decreasing in a (or p) and ^ rd is 
monotonically increasing in a (or p). Due to the opposite behaviors of the first- and second-hop 
SNRs as a function of a (or p), the optimal value of a (or p) minimizing the overall error 
probability is strictly between 0 and 1. 


D. Impact of the Modulation Alphabet Size M 
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In this subsection, we are interested in evaluating the performance of a variety of noncoherent 
signalings for M = 2,4, 8,16. To that end, we focus on a single-relay network with D 0r = 1.5 
(m) employing the near-optimum EH parameters, i.e., p = 0.8 for PS-NcDF and a = 0.4 for 
TS-NcDF. Note that p = 0.8 and a = 0.4 minimize, respectively, the SERs of PS-NcDF and 
TS-NcDF, as shown in Fig. |U However, it is difficult to analytically determine the optimum EH 
parameters. An interesting phenomenon observed is that while M-DPSK outperforms M-FSK 


for M < 4, the performance gain diminishes as M increases, as shown in Figs. |5(aj] and |5(b) 


In particular, for M — 4 shown in Fig. |5(b)[ the performance gain of DPSK over FSK is almost 
negligible. As M increases to M > 8, we see that M-FSK achieves a significant gain over 
M-DPSK. In particular, for a target SER of 10~ 2 , the performance gain of M-FSK over M- 


DPSK is approximately 4 dB and 8 dB for M = 8 and 16, respectively, as shown in Figs. 5(c) 


and |5(d)[ This implies that M-FSK achieves a promising energy gain over M-DPSK for large 
M. However, it is important to note that this energy gain comes at the cost of higher required 
bandwidth. In particular, the minimum double-sideband (DSB) bandwidth at an intermediate 
frequency (IF) assuming ideal rectangular Nyquist filtering If33l is W = 1/T S for M-DPSK and 
W = M/T s for M-FSK (M times that of M-DPSK). Nevertheless, considering the significant 
energy gain, M-FSK may be a more suitable solution, especially when M > 8, for EH relay 
systems due to the limited energy. 


VII. Conclusions 

We have proposed a noncoherent SWIPT framework embracing both PS-NcDF and TS- 
NcDF in a unified form, supporting M-FSK and M-DPSK. Following this framework, we 
developed exact noncoherent MFDs for PS-NcDF and TS-NcDF in a unified form, which 
involves integral evaluations. Nevertheless, the exact MFDs are useful for characterizing the 
optimum performance benchmark for noncoherent SWIPT. Furthermore, to make the MFDs 
suitable for practical implementation, we developed closed-form approximate MFDs achieving 
near-optimum performance. Numerical results demonstrated a performance tradeoff between the 
first and second hops through the adjustment of the time switching coefficient or power splitting 
factors. Moreover, the optimal values of these system parameters corresponding to the minimum 
SER are strictly between 0 and 1. Finally, we demonstrated that M-FSK is an energy-efficient 
solution for noncoherent SWIPT, especially when M > 8. 
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APPENDIX A: PROOF OF Lemma □ 

The conditional distribution of X 0 given X 2 = z is X 0 \x 2=z ~ CJ\f(0, S 0 ), where 

12 

,\ f, 1 I/, 

S 0 = 


+ Do, 

Di \z 2 c* 

| 2 c, 

Di \z 2 |c| 2 + D 0 


It can be shown that det(X 0 ) = D 0 Di|,£| 2 (1 + |c| 2 ) + Dq an d 




det(E 0 


Di|z| 2 |c| 2 + D 0 , —Di|z| 2 c* 

—Di|z| 2 c, Di|^| 2 + D 0 


After some manipulations, the conditional PDF is written as 


fx 0 \X 2 =z( x ) - 


/ I |2 \ pxn I_ np(i+lc| a ) _ 

exp P ( 


7T 2 Dq 


i + S(i + M 2 )M 2 


(A.l) 


(A.2) 


(A.3) 


The unconditional PDF is fx 0 ( x ) — J'Z=o f I, |2 _ (*)/|x 2 | 2 ( x )d 2 :, which is evaluated to 

X 0 | X 2 1 —-3? 


/*<>(*) = 


laia-CTir 
e ^oU+H 2 ) 


n 0 (i+|cp) 

1+ iilE2 (1+ | c | 2 ), 


vr 2 D 2 J 0 l + ^(l + | c | 2 )t 


e df. 


(A. 4) 


\xi+c*x 2 \ 2 

Applying change of variable x <- — 

1 ' f! 0 0+l c l ) t 


in (IA.4I) and after some manipulations, we have 


1^2 — cx i 


n 0 


]xi+c* x 2 \- 


fx 0 (x) = 


e n o( 1 +l c l 2 ) e nin 2 (H-|e|2) r f! 0 (i+|cp) 1 I x+ 


ir'-n$ agi(i + Id 2 ) J 0 


\x 1 + c*x 2 \ 


n 1 n 2 (i+\c\‘ 2 y 2 x J 


X 


(A.5) 


which is the same as (fT8l) upon some appropriate parameterizations. 

The derivation of the PDF of To follows similar lines as the derivation of fx 0 (x), which is 
briefly summarized as follows. The conditional PDF of Y fl given AT = x is given by 


M 


/ 


(y) = 


Y 0 \\X 2 \ 2 =x ya ' (ttD 0 ) M 




9=1 

<?Ap 


i + ii- 


exp — 


\y P 


Dq T Dix 


(A.6) 


Similarly, the unconditional PDF is obtained as 


fYoiy) = 


(ttd 0 ; 


M 


n«ir-i y P \ 2 

n 0 


i + i 22 


exp — 


\y P 


Dq T D^X 1 D 2 


1 __ 2 _ 

-e n 2 dx. 


(A.7) 


Using similar change of variable in (IA.7I) . we can obtain the final expression of (fT9l) . 
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APPENDIX B: Proof of Lemma [2] 
It follows that Y 2 ~ OV’fO. S 2 ), where 

1 + 7 , 7 c" 

= D n 


(B.l) 


7 c, l + 7|c| 2 

where 7 = Di/D 0 . Applying the property of the multivariate Gaussian distribution, the condi¬ 


tional distribution of Y 2 given Y\ follows a Gaussian distribution, i.e., Y 2 \ Yl = yi ~ CJ\f(tf 0 yi, a\ 
where 0 0 = jt~ and <r ( 2 = D 0 (l + yr-| c ) ■ Let Z 0 = Y 1 *Y 2 . It follows that Z 0 |y 1= . 


OJ’ 


CJ\f{do\yi\ 2 , al\yi\ 2 ). Thus, the unconditional PDF of X 0 is 


expf^ 

fz 0 (z) = _ 2 V g ° 


7T(To(Do + ^ 1 ) Jo X 


— exp 


Do T Di 


+ 


14 


CTn 


X 


°0 X J 


dx. 


According to ||29l eq. (3.471.9)], (IB.3I) is evaluated to 


/z » (2) = ^! +SJl) exp 


+ 


|4| 2 \z\ 


(B.2) 


(B.3) 


Dg + Di (Tq ctq / 

where Kq(x) is the zero-order modified Bessel function of the second ki nd ||29l eq. (8.432.1)]. 
Using the cartesian-to-polar transformation, /i' 0 e jO ° = Z 0 , the joint PDF of R 0 and 0 O is 

’2|44 


f Ro,@o (+0; 4) — 


7 T(Tq (Do + D] 


-r 0 exp 




cos(4 — Xc) 


Kn 


2 r 0 

(Tq V Dq + Di 
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14 


CTn 


(B.4) 


Since 0 = 0 O — Xc, the PDF of 0 is fe(9) = J°° /X 0 ,e 0 (r 0 , 6* + Xc)dr 0 , which evaluates to 

'2|4ko 


/e(4 = 


vra^Do + D x ) 7 0 


r 0 exp 


(Tn 


cos 9 I K 0 I —- 


( 7 q I/ Dq + Dx 


+ 


14 




It is verified by Mathematica® that the following equation hold+ 


xe ax K 0 (bx)dx = 


a arccos 


b 2 - a 2 ( b 2 - a 2 ) 3 / 2 


dr 0 . (B.5) 


(B.6) 


for b > 0 and St(a + b) > 0. Letting a = — 2 l'®j r ’° cos 9 and b = 2 


solve the integral of (IB .5b in closed-form as 
1 


fe(9) = 


2vr 


1 + gin 2 £ 


1 + 


cos 9 


Sin2fl + PM 


&0 V ^ 0+^1 


arccos — 


1 + (|b. 6I) . we 


COS# 


sin2g + \M*(t+Sh) - 


Finally, substituting $0 and a 2 } into the above expression yields the final expression in (CHT) . 


4 An alternative closed-form solution of the integral in dB~6t is available in f29l eq. (6.621.3)]. However, that expression is 
given in terms of the Gaussian hypergeometric function, which is more complex than the expression in ( fB~6t . 
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APPENDIX C: Proof of Lemma [4] 


From (1231) . we can rewrite I(e,/3) for (3e > 1 as 


/M) = ^°(Vf j + 


vT 


1 _(*+!) _± 
—e V e /e 


dx. 


ex 


(C.l) 


Over the integration interval (1/e, y//3fe) in (1C.II) . we have 0 < ^ < 1, 0 < e - ( x+ «) < 1, and 
0 < e - ^ < 1, meaning that the overall integrand is a finite bounded function. Moreover, the 
integrand is always between 0 and 1, which ensures fast convergence when approximated by 
the Gauss-Legendre quadrature (GLQ) lf34l . Note that the GLQ is one of the most efficient and 
accurate non-adaptive quadrature rules for a given number of arithmetic operations |[34l . Applying 
the JV-th order GLQ rule of 11301 eq. (25.4.29)] in (1C.lb yields I(e, /?) = /^(e, (3) + Rn, where 
In{£,/3) = + ^n(\, V^f)] an( ^ is given in (1251) . This gives the third case 

of (l24l) . For j3e < 1 in (1231) . it is easy to verify that the overall integrand is not necessarily smaller 
than one in this case. Thus, simply applying the GLQ rule for /3e < 1 in (1231) might give an 

approximation with slower convergence speed compared to the GLQ approximation for (3e > 1. 

1 1 

However, if we further enforce — < 1, then the pre-bracket coefficient — in (l23l) can scale 

down the approximation error for the integral, which can significantly accelerate the convergence 

1 

speed for the approximation. Thus, for (3e < 1 and — < 1, the N-th order GLQ rule of |[30l eq. 
(25.4.29)] readily applies, which yields I N (e,/3) = This is the 

same as the second case of (l24l) . For the remaining case with (3e < 1 and ^ > 1, we start with 
the alternative form of J(e, 3) in (l22l) . which can be written as J(e, 3) = E [ 1+ ^ Xo exp (— 1+ f Yo )] 
where the expectation is taken with respect to X 0 and X 0 is an exponential random variable with 
unit mean. Thus, the integral I(e,/3 ) can be interpreted as the statistical mean of the random 
variable g(X 0 ) = | ' V() exp ( — 1+ ^ eXo ) • Using the A r -th order estimate of the mean of random 
variables, we can approximate I{e,/3) by I N (e,(3) as follows If26l eq. (5.85)]: 
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In{£, ft) — ^ 


3 n & n g{x) 


n=0 


n\ 


dx 


(C.2) 


x=E{X 0 ]=1 

where /x rf = E[|X 0 — E\X () ]\ n ], n — 0, ■ ■ ■ , N. Using the Binomial theorem and the exponential 
distribution of X 0 , it is not hard to show that /i n = ^p-n!. The function g(x) can be 

represented as a composite function with g(x) = fo[go(x)}, where / 0 (f) = te~^ and t = g 0 (x) = 
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1/(1 + ex). Thus, Faa di Bruno’s formula If35l is applied as follows: 


d n fo[g 0 (x)\ 


dx 


n\ 


/ j n 

x=l nr=i irrii=n ]/[ 171 j 




n 


k =1 


a ( o k \x) 

k\ 


m k 


(C.3) 


X=1 


j =1 


where r = g 0 ( 1), = ]P” =1 rry, and Z + = {0,1,2, •••}. It can be shown that fo\t) = 

(q — /3t)(—/3y- 1 e-P t and g^\x) = ("+-^+1 • Applying these relationships in (IC.3I) . substituting 
the expression into (IC.2I) . and after tedious manipulations, we obtain the first case of (l24l) . 


APPENDIX D: Proof of Lemma [5] 

According to the MLD rule in (l27l) . the phase © r = , (l)yor(J) is a sufficient statistic for 

the MLD. Let V m = {e\6 m - < 9 < 9 m + denote the decision region for the message 

m, where 9 m = m = 0, • • • , M — 1. Then the ML decision rule is in favor of the message 
m r if O r e T> mr , m r — 0, ■ ■ ■ , M — 1. Thus, the transition probability is written as 


Pr(m r |m) = Pr 


[2(m r - m) - 1]tt < @ ^ < [2 (m r - m) + 1]tt 


m 


(D.l) 


M ~ ~ M 

The unified system model for T r in (U3al) can be rewritten as y 0r = h 0r [l,e jSrn ] T + n' 0r , 
where h 0r = cr 0r y / 7 (>^ors(/ — 1) and n' 0r = cr 0 r n 0r . It follows that h 0r ~ cx^y 0r ) and 


n or ~ 


CAf(0, <jQ r ). Applying the result in Lemma El the PDL of 0' r = 0 r — 9m is given by 


fe>(0) = 


27r 1 1 ~<ar sin2 9 

1 “ r l+2 70r 
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70 r cos 9 


\J 7or sin2 6 + 2 7or + 1 


arccos 


70 r COS 9 

70 r + 1 


(D.2) 


Thus, the transition probability in ( ID. II ) can be written as 


Pr(m r |m) = 


2(7717- —m) + l 
* M 7 


! 2( 


—m 


)-i 


2(771— m r ) + l 

r m 1 

r 2(771—771 t") — 1 __ 
M ^ 


fe'(9)d9 


f &r (a)da 


= Pr(m|m r ), 


(D.3a) 

(D.3b) 

(D.3c) 


where (ID.3bl) is due to the change of variable a i - 9 and /©/,(#) = /©^(—0). Note that the 

symmetry Pr(m r |m) = Pr(m|m r ) can be also attributed to the rotational invariance property of 
the PSK signal constellation. According to (ID.3I) . the transition probability can be calculated as 

21 7717—7711 -|-1 

f M 7r 

Pr(m r |m) = / f & (9)69. (D.4) 

/ 2|77lr — 7711 — 1 
J -M- 7r 

Finally, substituting (ID .21) into (ID.4I) yields the final expression in (|3TI) . 
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APPENDIX E: Proof of Lemma [6] 

Noticing that the transmitted PSK phase will most likely be mistaken for its two closest 
neighboring phases with equal probability PTll . we can approximate the transition probability as 
follows: for the two closest error transitions with \m r — mn\ = 1 or iff — 1, we have Pr(m r |m) ~ 
For the correct symbol transition with m r = m, we have Pr(m r |m) = 1 — SER. For all 
other error transitions, we simply ignore the probabilities. This approximation immediately gives 
(l32l) . It remains to derive the SER for M-DPSK. For M = 2, a closed-from exact expression of 
the SER is available, namely, SER = e? = - 7 -j-T— . For M > 4, however, no closed-form exact 
SER expression is reported in the literature for fading channels. To proceed, we resort to an upper 
bound on the SER for AWGN channels in If36l eq. (8.95)], which is the best known bound for 
AWGN channels. Taking the statistical average of this bound over the Rayleigh distribution by 
means of 11291 eq. (6.283.2)], we obtain a novel upper bound on the (average) SER of M-DPSK 
as £]? = 1.03^/^f [1 - . Overall, the SER of M-DPSK is given by 63- 


APPENDIX F: Proof of Theorem Q] 


Applying Total Probability Theorem and Lemma [5] the likelihood function (LF) conditioned 
on message m for the unified M-DPSK model in (fl3l) can be written as 

K M —1 

f {yod, {yrd}r = 1 \m) = f(yod\m) nz f {yrd\'W‘r)'P\m r —m\ (TOr) ■ (F.l) 

r =1 m r =0 

From (I13bl) . we see that y 0 d|m ~ CJ\f(0, £ orf ), where 
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After some manipulations, we have 
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From (I13cl) . we see that applying Lemma [H immediately yields the LF f(y r d\m r ) as follows: 
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Substituting (IF. 31) and (IF.4I) into (IF. II) and using the definitions of /3 0 (m), 6f(m r ), and /3~(m r ) 
given in Theorem [U we obtain the MLD in (f34l) for M-DPSK. For M-FSK in (fl4l) . following 
similar procedure as in (IF. 1 1) for M-DPSK, the LF can be written as 

K M—l r -i 


f(yod,{yrd}?=i\m) = f(yod\m) iie 


r=lm r =0 L 


(1 - e r )f(y rd \m r = m) + 
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, (F.5) 


Furthermore, the PDF of y 0 d in (I14bl) . conditioned on the knowledge of m, is 
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Applying Lemma Q] the PDF of y rd in (I14cl) given m r is given by 
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Finally, substituting (IF.61) and (IF. 71 ) into (IF.51) yields the MLD in ([331) for M-FSK. 
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